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Absolute proper motions are determined for stars and galaxies to V=17.5 
over a 450 square-degree area that encloses both Magellanic Clouds. The proper 
motions are based on photographic and CCD observations of the Yale/San Juan 
Southern Proper Motion program, which span over a baseline of 40 years. Multi- 
ple, local relative proper motion measures are combined in an overlap solution us- 
ing photometrically selected Galactic Disk stars to define a global relative system 
that is then transformed to absolute using external galaxies and Hipparcos stars 
to tie into the ICRS. The resulting catalog of 1.4 million objects is used to derive 
the mean absolute proper motions of the Large Magellanic Cloud and the Small 
Magellanic Cloud; (/x Q cos 5, ^s)lmc = (1.89, +0.39) ± (0.27,0.27) mas yr" 1 and 
(fi a cos 5, fis)sMC — (0.98,-1.01) ± (0.30,0.29) mas yr -1 . These mean motions 
are based on best-measured samples of 3822 LMC stars and 964 SMC stars. 
A dominant portion (0.25 mas yr^ 1 ) of the formal errors is due to the esti- 
mated uncertainty in the inertial system of the Hipparcos Catalog stars used to 
anchor the bright end of our proper motion measures. A more precise deter- 
mination can be made for the proper motion of the SMC relative to the LMC; 
{Ha COS 5,Hs)sMC-LMC = (-0.91, -1.49) ± (0.16,0.15) mas yr -1 . This differential 
value is combined with measurements of the proper motion of the LMC taken 
from the literature to produce new absolute proper-motion determinations for 
the SMC, as well as an estimate of the total velocity difference of the two clouds 
to within ±54 kms -1 . The absolute proper motion results are consistent with 
the Clouds' orbits being marginally bound to the Milky Way, albeit on an elon- 
gated orbit. The inferred relative velocity between the Clouds places them near 
their binding energy limit and, thus, no definitive conclusion can be made as to 
whether or not the Clouds are bound to one another. 

Subject headings: astrometry — catalogs — proper motions — Magellanic Clouds 



1. Introduction 

The Magellanic Clouds have provided astronomers with a wide variety of information, 
from the small to the large hierarchy of objects in the Universe. They are the first step in the 
cosmic-distance ladder, as well as a proxy for the low-metal gas-rich galaxies assembled in 
the early universe, they have zones of recent massive star formation, young 1-3 Gyr globular 
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clusters and a variety of old pulsating stars and planetary nebulae. They contain substantial 
amounts of gas and dust in a violent interstellar medium, harbor the closest supernova in 
recent years, and are being used as a testbed for dark matter searches. 

The Magellanic Clouds are also the prime example of a galaxy-galaxy interaction, based 
on several lines of evidence, close e nough to be studied in detail: An ap parently starless 
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Arm precedes the clouds (IPutman et al.lll998l ); an HI envelope surrounds the whole system 
and a collection of High Velocity Clouds seems to be "raining" over the Galaxy ( IOranoll2004l ). 
Despite such convenient observational circumstances, on the theoretical side no dynamical 
model has been able to reproduce all these phenomena simultaneously. 

It is widely accepted that the complexity and intricacies of the Magellanic Clouds' 
external and internal features, have been largely determined by the orbit they have followed 
in the past few Gigayears. Due to their large distance, about 50 and 60 kpc to the LMC and 
the SMC, respectively, only radial velocities were precise enough to provide some assessment 
of their spatial velo city In fact, line-of-sight measurements of the Clouds began about a 
hundred years ago (I Wilson! Il 9151 ). while proper motion measurements of a useful accuracy 
were possible only in the 1990's. 
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in which the Magellanic Clouds were bound to each other and to the Milky Way. Such 
scenario relied heavily on the fact that the Galactic gravitational potential used (isothermal 
sphere) yields such results by default, and proper motion errors were not small enough to 
refine the tangential velocities. 
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More interestingly, HS T-based results kallivavalil et al.lbood . l2006bl Jal: IPiatek et al1l2008 



Kallivavalil et al.l 120091 ) coupled to more modern and cosmologically inspired dark matter 
Halo models, suggest that the Clouds are not bound to the Galaxy, opposite to the long-held 
paradigm. 



Twenty years have passed since the very first proper motion measurements of the Clouds, 
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and it is only now with the Yale-San Juan Southern Proper Motion (SPM) program - briefly 
explained in Section [2]- that for the first time a wide- field astrometric proper motion survey 
of the Magellanic System is finally completed. All the lst-epoch (early 1970's) and part of 
the 2nd-epoch (early 1990's) SPM material used in this work are photographic plates. Their 
processing is briefly summarised in Section [3j but a more detailed explanation can be found 
in Girard et al. (2010). 

A substantial part of our 2nd-epoch data comes from SPM CCD observations. A short 
explanation of the data acquisition, quality control and processing is explained in Section HI 
To supplement our 2nd-epoch plate data with CCD-quality positions, we haye inc luded mean 



positions at Julian epoch 2000.0 from the UCAC2 catalog (IZacharias et al.ll2004[ ). Section [5] 



explains how these data were selected and included in this investigation. 

In contrast to other SPM reductions, as explained in Section El relative proper motions 
measured in CCD-size fields of view were combined into a single common extended and 
accurately defined global system. Although our zero point accuracy, i.e. how well our 
reference frame is linked to the International Celestial Reference System (ICRS), is ultimately 
limited by Hipparcos accuracy, our very precise relative proper motions over the whole field 
of view, enabled us to measure the proper motion of the SMC with respect to the LMC, at 
a precision comparable to the quoted errors of space-based proper motions. 

It is this capability that we exploit to obtain new measurements of the proper motion 
of the SMC based on previously published LMC proper motions. Section [7] contains the 
main results of this paper regarding the proper motion of the Clouds, absolute and relative. 
Section [8] has a discussion of the implications of our results on the current understanding 
of the dynamics of the Magellanic System. Finally, Section [9] states the conclusions of this 
paper, and future plans already in consideration to improve the current results. 



2. The SPM program 

This investigation is part of the SPM program, a joint venture of Yale University and 
Universidad Nacional de San Juan in Argentina. The SPM program was initiated in the early 
960s by D. Brow er and J. Schilt as a joint enterprise of the Yale and Columbia Universities 



(jWesselinkl 1 19 741 ) . The goal of the SPM program is to provide absolute proper motions, 



positions, and BV photometry for the Southern sky to a limiting magnitude of V ~ 18. 

The SPM program makes use of the Yale Southern Observatory's double astrograph at 
Cesco Observatory in the foothills of the Andes mountains in El Leoncito, Argentina. This 
telescope consists of two 51-cm refractors, designed for photography in the blue and yellow 
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bands, respectively. The first-epoch survey, taken between 1965 and 1974 was made on glass 
photographic plates, exposed simultaneously in blue- yellow pairs and always centered on the 
meridian. The plates' field of view (FOV) extends over an area of 6.3° x 6.3°. The sky south 
of 5 = —17° was observed in the first-epoch period. 

Second-epoch SPM plate observations were begun in 1988. By the mid 1990's, with 
only a third of the second-epoch survey completed, Kodak discontinued the production of 
the photographic plates. In 2000, with funding from the NSF, a CCD camera system was 
installed on the double astrograph to replace the photographic plate holders. A PixelVision 
4Kx4K CCD camera (0.94° x 0.94° FOV) was placed in the yellow lens focal plane, and an 
Apogee Ap-8 lKxlK (0.37° x 0.37° FOV) was fitted in the blue focal plane. In 2004, the 
Apogee IK camera was upgraded to an Apogee Alta E42 2Kx2K (0.42° x 0.42° FOV) with 
funds from the Argentine CONICET. In order to achieve a limiting magnitude similar to the 
first-epoch plate material, the CCD survey consists of 2-minute exposures in both cameras. 

In the past decade, catalogs of the SPM program covering various parts of the sky 
have been published, as th e second-epoch ma terial bec ame available for its as trometric re- 
duction. Catalogs SPM1 Jpiatais et alJll998h . SPM2 (Ivan Altena et al.lll999h . and SPM3 
(IGirard et al.l 12004 ) are based on photographic plates only, in both first and second epochs. 
The plates were scanned either with the Yale PDS (using input lists of selected objects) or 
the USNO Precision Measuring Machine (PMM) for the whole plate. Also, two different 
centering algorithms have been used to meas ure the image centers in the sc ans, the Yale 



2D-G aussian fit from iLee &: van Altenal (119831 ). or the USNO circular fit from lMonet et al. 
(l2003h . 



Since 2004, regular CCD observations have been carried out to finish the second-epoch 
survey of the SPM program. By December 2008, the survey was effectively completed for 
the sky south of 5 = —20°. Subsequently, the SPM4 catalog, based on all available plate 
and CCD data, was completed in late 2009 and is currently available, (Girard et al. 2010). 
SPM4 includes ~100 million objects south of 5 = —20°, brighter than V ~ 18. Many of the 
data-processing procedures, software and protocols developed in this investigation were also 
used in the construction of the SPM4 catalog. 



3. The Plate Data 



3.1. Observations 



The SPM survey fields are on 5° centers in declination and a maximum of 5° separation 
in right ascension, providing at least a full degree of overlap between adjacent plates. Each 
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17-inch x 17-inch plate covers an area of 6.3° x 6.3° (55.1 "/mm plate scale) and consists 
of a 2-hour and a 2-minute offset exposure. All observations were made with a wire grat- 
ing over the objectives, producing measurable diffraction images out to third order. The 
grating constant is 3.8 magnitudes, thus, along with the offset short exposure, effectively 
increasing the dynamic range of each plate allowing measurement of external galaxies and 
bright Hipparcos stars. Fields were observed simultaneously in blue and yellow passbands, 
and there are some fields with repeated blue and/or yellow plates from the same epoch. See 



Girard et al.l ( 120041 ) for a more complete description of the SPM plate material. 



Tabled] lists the SPM plates used in this investigation. During the course of this research, 
it was found that some plates yielded unusually deviant results. A visual examination of 
the suspect plates revealed that the stellar images suffered from significant defects, possibly 
caused by poor guiding, polar misalignment, or poor focus. These plates, 0750B, 0751B, 
1371B, 1357Y and 1373Y, were therefore discarded for the research presented here. Coverage 
in these areas was not affected, since only one of the two plates at a given epoch per field 
affected was discarded. Figure [1] shows the distribution on the sky of the plates used in this 
work. Twenty two (22) SPM regions were studied, of which seven (7) have 2nd-epoch plates. 
For the areas with 2nd-epoch plate data, no 2nd-epoch CCD observations were made. 



3.2. Astrometric Reduction and Photometric Calibration 



The SPM plates used in this investigation were scanned with the Precision Measuring 
Machine (PMM) o f the USNO Flagst aff Station. For more details about the PMM setup 
and operations see iMonet et al.l (120031 ) . In a collaboration betwe en the USNO- Washin gton 
and the Yale Astrometry Group, the StarScan reduction pipeline (jZacharias et al.ll2008l ) was 
modified to analyze the PMM pixel data of the SPM plates to produce a list of detections, 
image centers and photometric indices. The astrometric and photometric reductions then 
proceeded as follows: 1) cross-identification of detections to an input catalog, including 2) 
recognition and identification of central and higher grating orders; 3) photometric calibra- 
tion to obtain BV; 4) correction for Atmospheric Refraction; 5) correction for Magnitude 
Equation, which also combines grating-order systems; 6) transformation of short-exposure 
positions into the long-exposure system; and 7) astrometric solution into Tycho-2 to obtain 
(a, 5). 

The in put catalog referenced in step 1) is a c ompilation of the following e xternal catalogs: 



Hipp arcos (IPerrvman fc ESAl Il997f ). Tycho-2 flHog et all l2000h. UC AC2 (IZacharias et al. 



2004 ). 2MASS point-so urce and extended-so urce (ISkrutskie et al.ll2006), LEDA galaxies with 



DENIS measurements (IPaturel et ajj|2005j ), and the QSO catalog of 



Veron-Cetty fc Veron 
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( 120061 ). In order for an object to be included in our study, it must appear in one or more 
of these listed catalogs. A thorough explanation of all these procedures can be found in 
Girard et al. (2010). After the above processing, one has positions (x,y) properly calibrated 
into a common system within each plate, with computed positional errors and astronomical 
coordinates (a, 5) on the ICRS, as realized by Tycho-2. 



3.3. Evaluation of the plate data 

Well measured stars on the plates in all of the orders have positional errors between 
0.9 /zm and 1.6 fim (50 mas and 90 mas). These errors, which only assess measurement 
uncertainties, are consistent with the precision expected for a good centering procedure, 
based on previous experience with the SPM plates. 

A single final position (x, y) per star per plate is obtained, from the positional-error- 
weighted-average of the available measurements. As expected, the error of the final position 
and magnitude will depend on the number of grating orders contributing to their calculation. 
If average measurements from different plates were later averaged to obtain a final number, 
then other errors, random and systematic, would come into play, and the error budget of 
this final result should include the additional sources of uncertainty. 

As with earlier SPM catalogs, an approximate estimate of the relative completeness 
between the plate data and 2MASS, can be made. Figure [2] shows the V magnitude distri- 
bution of all stars detected on a yellow plate, compared to the Vjk magnitude distribution 
of all 2MASS point source stars in the same field. Vjk = J + 2.79( J — K) is an approxi mate 



empirical transformation from 2MASS JHK to V determined by iGirard et al.l (120041 ) and 
found to be valid for a relatively wide range of spectral types. It can be seen that the SPM 
plates have a completeness similar to that of 2MASS up to V = 17.5 and a falloff after that. 

In general, compared to previous SPM processing, these plate data have significantly 
fewer false detections, and a better correction of systematics in the detected positions asso- 
ciated with the scanning process. 



3.4. SPM 2nd-epoch plate data 

As explained in subsection [371], seven SPM fields have 2nd-epoch plates. For this reason, 
no 2nd-epoch CCD observations were made at these locations. In order to facilitate the 
managing of files and software, the 2nd-epoch plates were divided in CCD-size fields, to 
emulate the overlap scheme of the 2nd-epoch CCD frames (See Section [472]) . This way, all 
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the 2nd-epoch material, regardless of its type, would have a uniform format and structure, 
for programming purposes. From the about 90 CCD pointings that usually cover one SPM 
field, half were used to divide the blue plate and the other half to divide the yellow plate, 
in such way that these yellow and blue fields overlap in a similar way as do the real CCD 
frames. 



4. The CCD data 

4.1. CCD cameras 

The CCD system consists of four cameras, one each for the blue- and yellow-optimized 
lenses and two focus-sensor cameras, again, one for each lens. The CCD camera for the 
Yellow telescope is a PixelVision camera with a Cryotiger Chiller cooled (-85°C) 4K x 4K 
Loral chip with 15 \im pixels, which translates into a pixel size of 0."83 and a total area 
of 0°.94 x 0°.94 degrees. The pixel size is well matched to the YSO site where the seeing 
conditions usually yield an image FWHM of 2-3", corresponding to 3-4 pixels per FWHM. 
This is optimal sampling for the derivation of astrometric image centers, based on SPM 
experience with digital image centering. The unthinned and front-illuminated Loral chip is 
fitted with a fixed Custom Scientific V-band filter. 

The Blue telescope was first fitted with an Apogee AP-8 camera that utilized a thermo- 
electrically-cooled (-40°C) and back-illuminated IK x IK Site chip with 24 /j,m pixels, which 
translated into a pixel size of l."32 and a total area of 22. '57 x 22. '57. In May 2005, an 
upgrade was made by replacing this camera with a new Apogee Alta E42 back-illuminated 
2K x 2K chip, with 13.5 \im pixels, which correspond to a pixel size of 0."74 and a total area 
of 24. '8 x 24. '8. Centered on the same field as the larger PV yellow camera, the purpose of 
the blue CCD camera is to provide B-band CCD photometry for the stars that fall into its 
FOV. 

Data from the yellow PV camera were used for both astrometry and photometry, while 
the blue Apogee and Alta cameras observations were used only for the photometric calibra- 
tion of the blue plates. 

4.2. Observations 

As a norm for the SPM, CCD observations are done always within \ h ?>Q m of the meridian, 
in 2-minute exposures, with the wire grating placed so that the diffraction pattern is at 
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about 45° from the E-W line. Normally, an E-W orientation is ideal to avoid differential 
color refraction effects within the diffraction pattern, but in this case, a diagonal orientation 
prevents the saturated central-order image of a bright star from spoiling its grating images 
by either row or column bleeding. The CCD pointings conform to a two-fold overlap coverage 
scheme for the PV camera, as shown in Figure [TJ About 90 PV CCD frames cover one single 
SPM field. Only targets that did not have 2nd-epoch plate data were observed with the 
CCD, except for a few special targets. 

A total of 1310 CCD pointings were observed for this investigation. Each pointing was 
planned to be observed only once, except for 90 of them extending over a 6° x 6° field around 
(a, 5) = (3 /l 44 m 33 s , -71°40'18"), within the area delimited by the bold black line in Figured] 
This area corresponds to the non-SPM VMC field, for the Variable Stars in the Magellanic 
Clouds study done from 1965 to 1968 by A. J. Wesselink, and comprises seventy 60-minute 
exposure blue plates without the wire objective grating, reaching a limiting magnitude of 
about 18. Numerous repeated CCD observations were performed on this area to provide 
suitable 2nd-epoch observations for this material. The VMC plates were not used for this 
work because they require their own special reduction, different from the one used in normal 
SPM plates, but the repeated CCD observations were indeed used for this investigation. The 
VMC CCD targets have at least 13 good observations each, with some of them having up 
to 20 good observations. 

The criteria used to qualify an observed frame as acceptable are: FWHM < 3." 5 , 
limiting magnitude Vu m > 17 and standard error < 120 mas from an astrometric solution 
into the UCAC2. If a frame fails any of these limits then it was taken again until it passed 
all of them. Nonetheless, all frames regardless of quality are saved and processed, and only 
later in the astrometric reduction are discarded, if they prove to be too bad for any use. In 
total, 5422 CCD frames were processed for this investigation. 



4.3. Astrometric Reduction and Photometric Calibration 

Data from all the CCD cameras went through the usual processing to calibrate the flux 
detected by the electronics, for the zero charge of the chip (bias), accumulated signal from 
the electronics dark current (dark) and different response to light from each pixel (flat). 
Details of these procedures can be found in Girard et al. (2010). 

Image detection on the processed CCD frames is done using SExtractor Version 2.4.4 



( IBertin & Arnoutslll996f ). from which a preliminary centroid and an aperture instrumental 



magnitude are read. SExtractor centroids are then used as input positions to compute more 



- 10 - 



precise centers, based on the Yale 2D gaussian centering algorithm (ILee fc van Altenalll983l ). 
Typically more than 90% of SExtractor detections are centered. A significant reduction in 
the rms of the positions for repeated observations of the PV frames was found when using 
the Yale-based centers as compared to the SExtractor-based centroids. 

The process of transforming the Yale-based centers (x, y) and SExtractor-based instru- 
mental magnitudes m inst into calibrated (a, 8) and BV, for the CCD frames, follows a 
procedure similar to that of the plates (Section 13.21) . with the exception that the external 
catalog used for the astrometric reduction of the CCDs is UCAC2 (Girard et al. 2010). As 
a result of these procedures, all the detected positions (x, y) are properly calibrated within 
each CCD frame, with computed positional errors and equatorial coordinates (a, 5) in the 
ICRS, as realized by UCAC2. 



4.4. Evaluation of the CCD data 

The single-image centering precision for well measured stars (V < 15) in the CCD 
frames is 0.5 jum (25 mas), worsening for the faintest stars where it reaches about 2 /im (100 
mas). A single final position (x,y) and magnitude per star per CCD is obtained, from the 
positional-error-weighted-average of the available measurements. The final positional and 
photometric errors will depend on the various image orders contributing to the final value. 

At this point, each star in each CCD frame has a master ID identification, a position 
(x, y) with errors, a calibrated B or V magnitude and a UCAC2-based (a, 5). Because stars 
were identified in the CCD frames in the same way as in the plates, similar charateristics 
regarding completeness were expected when compared to 2MASS (See Figure [2]). The CCD 
data show in general a completeness magnitude of about V = 18. 



5. UCAC2 CCD-positions as supplement for the 2nd-epoch plate data 

As seen in Figure [IJ some fields on and adjacent to the Magellanic Clouds have only 
photographic plates as 2nd-epoch material. Given the lower quality of plate images compared 
to CCD observations, the measured proper motions in these plate-only fields, will have 
significantly larger errors than those coming from the combination of plate and CCD data. 



In an attempt to counter this and achieve a more homogeneous quality in the final 
measurements, we decided to supple ment the 2nd-epoch p late measures with epoch 2000 
positions from the UCAC2 Catalog ( iZacharias et al.l 120041 ) . These are mostly based on 
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CCD observational with the USNO 8 inch (0.2 m) Twin Astrograph from Cerro Tololo 
International Observatory in Chile. The UCAC2 data were collected in such way that it 
mimics our 2nd-epoch fields. The precision of the positions are 15 to 70 mas, depending on 
magnitude, with claimed estimated systematic errors of 10 mas or below. UCAC2 provides 
only crude magnitudes in a single nonstandard bandpass R\jcaci between V and R, and its 
limiting magnitude is about R ~ 16. 

A significant number of faint stars in our 2nd-epoch SPM plates do not have a counter- 
part in the UCAC2 catalog. This means we are only sampling stars about V < 16 for these 
SPM fields when using UCAC2 positions. UCAC2 completeness compared to 2MASS and 
SPM data, can be seen in Figure [2j 



6. Obtaining the proper motions 

Given that a substantial part of the 2nd-epoch SPM material used in this work comes 
from CCD frames, previous SPM procedures for obtaining proper motions (used on 1st- 
and 2nd-epoch plates), could not be straightforwardly applied. Since the CCD's FOV is 
about 40 times smaller than the plate's FOV, the number of reference stars available to 
measure relative proper motions in each field is proportionally smaller. A simple cut in 
magnitude, as in past SPM reductions, would result in too few reference objects per CCD 
frame; in particular, the extragalactic objects needed to transform relative proper motions 
into absolute ones. 

However, if we select reference stars belonging to some specific population of the Galaxy, 
it is reasonable to assume that they have a mean absolute motion along an extended area 
on the sk}|f| that can be parametrized as a smooth function of (a, 5). Moreover, within a 
CCD FOV, their mean motion has a very small gradient, if any. It then becomes a matter 
of precisely quantifying, over the whole field of view, the measured mean relative proper 
motion of all known extragalactic objects, which is simply the reflex of the mean absolute 
proper motion of the reference stars. Applying such a function to the measured relative 



2 Strictly speaking, the supplemental UCAC2 epoch-2000 positions we use include up to a few years' 
worth of UCAC2 proper motions. For the areas of sky in this study, the UCAC2 CCD observations were 
made around 1998. Thus the positions we employ as supplements are practically those of the UCAC CCD 
program. 

3 Stars behave like a collisionless system, all moving under the influence of the same general background 
gravitational field, therefore we can expect them to have a global smooth distribution in their velocities, 
with some scatter around a mean value at any given location. 
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proper motions converts them to absolute. 

The layout of the CCD fields, where substantial overlap exists (20% to 50%), contains 
a wealth of information that can provide linkage of the reference system across the observed 
region of the Clouds that is limited only by measurement errors. A key point of this investi- 
gation is therefore to find a procedure to utilize this large overlap to produce a precise global 
reference system. 



6.1. Relative proper motions in the small CCD FOVs 

After all the plate and CCD data processing described above, the following data are 
ready for the measurement of proper motions: 

• lst-epoch positions on the plates, which have been corrected for systematics, as much 
as Tycho-2 precision and number of stars available allows. 

• 2nd-epoch positions on SPM CCD frames, SPM plates, or from UCAC2 . CCD po- 
sitions are mostly free of systematic errors. Plate positions have been corrected for 
systematics, as much as Tycho-2 stars precision and number of stars allows. 

Henceforth, we will refer to each of the 2nd-epoch CCD-size frames as a brick, regardless 
of the source of its data (SPM plate, SPM CCD or UCAC2). For each 2nd-epoch brick, the 
corresponding lst-epoch plate area was reprojected onto a tangential plane centered on the 
brick. Then a quadratic solution was computed to transform the reprojected lst-epoch 
plate's (x, y) into the 2nd-epoch brick's (x, y), forcing the chosen reference stars to have zero 
mean proper motion. The quadratic terms in the solution were meant to model systematic 
errors, either from the plate (uncorrected Optical Field Angle Distortion) or from Tycho-2 
proper motion systematic^ unavoidably propagated backwards in time into the computed 
lst-epoch (x,y). These solutions yield measured relative proper motions in each brick. The 
vast majority of the solutions only needed linear terms and the typical standard error of 
the solutions varied from about 5 to 8 mas yr _1 , which is dominated by the intrinsic proper 
motion dispersion of the reference stars. 

The reference stars were chosen with the followi ng criteria: 1 < V — J < 1.5 , 0.25 < 



J — K s < 0.65 and 13 < K s < 15, which according to iNikolaev fc Weinberg! (120001 ) isolates 



4 Such systematics were indeed found later on in some tests. The Tycho-2 Catalog, although an astrometric 
catalog based on late epoch space-based data, has early epoch positions from ground-based data that are 
known to suffer from significant magnitude equation. 
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mostly G-M dwarfs in the Galaxy disk, located between 0.4 to 1.6 kpc from us, with an 
estimated mean distance of ~ 650 pc. Their distribution in V magnitude ranges mostly 
from V — 15 to V — 16.5. In general there are between 200 and 500 reference stars per 
field, depending on the galactic latitude (our fields extend from b ~ —50° to b ~ —20°). The 
intrinsic proper motion dispersion of the Galactic disk reference stars was seen to increase 
with |6|, reflecting the changing kinematics of the Galaxy along it. 

Despite the photometric cuts to select the reference stars, contamination by LMC and 
SMC stars could not be avoided in the densest parts of the Clouds. This forced us for the 
time being to restrict our investigation to those fields in which we trust the relative proper 
motions, as being measured with respect to bona fide Galactic foreground stars. These areas 
were defined as shown in Figure [31 Although a substantial number of Magellanic Clouds 
stars were lost in this selection, on the other hand, confusion due to image crowding at these 
locations render these fields useless anyway, due to the risk of misidentifications. 

From the initial 13880 bricks available, 12180 are in the non-contaminated areas. After 
rejecting bricks from the discarded plates (Section l3.ip . 10900 are left to build the catalog 
of proper motions, mostly outside the Magellanic Clouds. In order to increase the number 
of LMC and SMC stars measured at the end, contaminated fields that overlap with this 
non-contaminated catalog were later on directly tied into it, and common stars had their 
proper motions averaged. 

6.2. Combining the proper motions 

Once relative proper motions with respect to the Galactic foreground stars had been 
measured in the non-contaminated area, different approaches were tried to combine them 
into a single global well-defined reference frame. Reference stars were chosen hoping that 
their mean motion along the sky could be described by a smooth function. This goal was 
indeed attained as confirmed by the fact that the measured relative proper motion of all 
known extragalactic objects were very precisely fit by a quadratic polynomial in (a, S). 

At this point, applying this polynomial to the relative proper motions to convert them 
to absolute, and then averaging all measurements, is a way to combine all the information 
available per star. But this would yield a catalog with a rather noisy zero point as one moves 
along the sky, mostly due to the real intrinsic proper motion dispersion of the reference stars 
and the fact that two overlapping frames may not have the same reference stars. Given two 
frames with about 50% overlap between them, means that both have 50% of the reference 
stars in common, while the other 50% are different. Given two samples with N data points 
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each, both with the same dispersion o and both having 50% of their points in common, it 

can be easily shown that their individual mean values typically differ by y~jjO~. In other 
words, we can expect that two overlapping frames typically differ in their relative proper 
motion systems by ~ 0.6 mas yr _1 , for N = 220 and o = 6.5 mas yr _1 . This is too large for 
the level of precision that we want to obtain. 



Eichhorn 


1960; 


Jeffervs 


1979 



Stockl Il98ll : lEichhornlll988l ; iTaffl 11988) was considered to link every frame's relative system 
into a global one. Each frame, containing a sample of the whole population of reference 
stars, realizes a local system that could deviate from the global system, due to statistical, 
measurement and/or systematic errors. A linear function per frame was considered sufficient 
to describe the difference between the frame and the global system. Therefore three para- 
meters per frame, per proper-motion component (fi a cos 5 and fis are solved separately) need 
to be determined. The coefficients of the quadratic polynomial that globally describes the 
mean motion of the reference stars are determined as well. To simultaneously solve several 
thousands frames, we would need to invert an approximately 40, 000 x 40, 000 matrix to 
get the parameters values and their errors. For this project, such a scheme was deemed 
impractical at this time. Nonetheless, ideas about using this approach over smaller areas 
first, like the SPM fields, and then performing another block adjustment solution to join 
these regions, are being adopted for future work. 

A more practical approach was considered to bring each individual frame's reference 
system closer to the global one. A single frame's reference system is statistically more 
deviant from the global system than is the system defined by all stars, since the proper- 
motion reference stars are only a subset of all stars. Using all stars in common between 
overlapping frames, we can adjust each frame's proper motions to agree with the average of 
the surrounding fields. This adjustment also helps to correct residual distortions, as they 
are statistically smoothed out in the average frame. To avoid frames drifting away from the 
global system as they are being aligned to one another, all reference stars in the field are 
explicitly assigned a relative proper motion of zero. Once the adjustments are applied, new 
averages can be computed, and the whole process is iterated until the adjustments converge 
to zero. 

Before making these adjustments, we first checked to see if systematic differences existed 
as a function of the magnitude. A non-negligible number of frames exhibited systematic 
trends with magnitude when compared to the average frame (See Figure H]). In general, a 
linear function of the coordinates in the field would take care of the geometrical distortion, 
but the magnitude equation required a smoothed-localized median, which can trace the 
general trend better than any parametrized fit. Therefore, each frame is first corrected for 
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its differential magnitude equation with respect to the average frame, and then we proceed to 
correct for the distortion, following the iterative procedure explained above. Each frame has 
typically about 1500 to 2000 stars, with some of them having up to 5000 stars, to compute 
the adjustments. 

With these improved relative proper motions, the quadratic polynomial that describes 
the mean reflex proper motion of the extragalactic objects was computed and used to trans- 
form the relative proper motions into absolute ones. Polynomials for fi a cos 5 and as 
functions of a and 5 were calculated separately, using proper-motion-error-weighted least 
squares. A total of 5351 external galaxies were used across the 450 sq-degree area and the 
formal errors of the polynomials, computed at the center of LMC and SMC, based on the 
full covariance matrix, amount to 0.03 mas yr^ 1 and 0.06 mas yr -1 respectively. 

At this stage, a final absolute proper motion and its corresponding error are obtained per 
star, from the error-weighted average of all the individual absolute proper motions obtained 
for it. A total of 1337050 objects in the non-contaminated fields had final proper motions, 
including a good number of LMC and SMC stars located in the outskirts of these galaxies. 
For such reason, it is named the outside catalog. 

In order to increase the number of LMC and SMC stars with measured absolute proper 
motions, contaminated fields that overlapped with the outside catalog just obtained, were 
directly tied into it by computing a linear solution to correct their relative proper motions 
and put them into its absolute system. A total of 678 additional fields were added with this 
procedure, eventually increasing the number of LMC and SMC stars by about 30% and 50%, 
respectively. 

6.3. Zero point global correction of the Absolute Proper Motions 

Since we are now theoretically on the absolute reference frame defined by the external 
galaxies, our catalog should be -within measurement errors- in the same reference frame 
system of other known catalogs of absolute proper motion. When we checked our measures 
(fi a cos 5, us) of 1356 Hipparcos stars against the Hipparcos Catalogue, we found a significant 
difference between the two, that amounts to 



t^ct cos S,Hipparcos f^a cos S,This work 0.49 i 0.07 maS yr 

Hs,Hipparcos — fJ>s,This work = —1-21 ± 0.07 mas yr 



-1 



(1) 

(2) 



-1 



The source of this systematic difference could be indeed in any (or both) of the two 
catalogs. Although Hipparcos is the most accurate optical astrometric catalog published 
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so far, it has significant correlations between the astrometric parameters (position, proper 
motion and parallax) of different stars, when they are less than about 5 degrees apart on 
the sky, and also between the astrom etric parameters for a g iven star, due to the special 



me asurement pr i nciple of Hipparcos (jPerryman fc ESAl 1 199 71 ). In fact, a new reduction 



by Ivan Leeuwenl (120071 ) was performed to correct some systematic correlations in the data. 
However no significant difference was found between the old reference frame (used in this 
investigation) and the new one. 

On the other hand, our plate measurements of the galaxies are not error-free. Tests were 
run to compare the final absolute proper motions if only galaxies with 2nd-epoch CCD data, 
or all of them, were used to compute the quadratic polynomials to transform relative proper 
motions into absolute. The polynomial for fi a cos 5 exhibited noticeable differences around 
the SMC. This was not completely unexpected given that at that location, four SPM fields 
have plate-only data, but in any case it points to the fact that galaxies in the 2nd-epoch 
plates may introduce problems. Hence, galaxies in the lst-epoch plates may do the same. 



As determined by iGirard et al.l (119981 ). the SPM plate material exhibits magnitude 
equation for galaxies and for stars that differ in functional form. However, the two could be 
brought into approximate agreement by adding an offset of -0.7 to the magnitudes of galaxy 
images before calculating the magnitude-equation correction as determined from stellar im- 
ages. In the present study, we have applied the same -0.7 mag offset to galaxy images for 
the purpose of magnitude-equation correction only. 

A comparison of the V magnitudes for the galaxies in our work, showed that plate pho- 
tometry returned signficantly fainter magnitudes than in CCD calibrations. Not suprisingly, 
this may have affected the magnitude equation correction in the plates, since it was the 
calibrated plate photographic magnitude that was used for such purpose. More surprising 
though, was to find that Eq. [TJ showed a linear trend versus a with a slope of about 15 
yuas yr _1 per degree (a bit smaller when using only galaxies with CCD data). The LMC 
and SMC center of mass positions are separated by about 20°, for which the trend above 
indicates a zero point shift of 0.3 mas yr _1 , that could in principle be related to Hipparcos' 
systematics. On the other hand, systematics in Eq. [2] versus 5 were also seen, that look to 
be related to the plates location and layout. 



The investigation of the SPM plates magnitude equation done by IGirard et al.l (119981 ) 
found that magnitude equation terms^l varied more or less uniformly from 1st to 2nd-epoch 
plates, so we can expect more or less uniform offsets with Hipparcos proper motions, if any 
residual uncorrected magnitude equation were still present in the data. Consequently, Eqs. 



A polynomial in (X, Y, m) is used to describe the magnitude equation correction. 
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[TJ and [2] were applied to all the absolute proper motions, putting our catalog on the system of 
the ICRS via Hipparcos. The inaccuracy of its zero point is dominated mostly by Hipparcos' 
systematic error of 0.25 mas yr -1 , since the quadratic polynomial, as defined by the external 
galaxies, was in general very accurate, being within 0.1 mas yr -1 of error for most of the 450 
square-degree area studied. 

6.4. Final Catalogue of Proper Motions - Evaluation of Errors 

The final catalog of absolute proper motions at this point has 1,448,438 objects, with 
the following data listed: a, 5, V; V — J, J — K s , H — K s when available; absolute proper 
motions \i a cos 5, [is and their formal errors e^ aCOS s,e^ s in mas yr -1 , number of data points 
used, number of data points rejected (outliers were rejected based on their normalized errors), 
a flag to indicate Hipparcos, Tycho-2, 2MASS extended sources, confirmed LEDA Galaxies 
and QSOs, a flag to indicate if the object is or is not a reference star, and the 2mass ID 
when available. The overall distribution of the stars in the catalog can be seen in Figure [5J 

In the final catalog of proper motions, stars with V < 12 have formal proper motion 
errors of about 0.5 mas yr -1 , and well measured stars with 12 < V < 15.5, have values 
that range from 0.5 to 1.3 mas yr -1 . These only reflect measurement errors, as they are 
based on the positional errors and the epoch difference. When proper motions for the same 
star measured in different bricks are averaged together, statistical and systematic deviations 
between the bricks must be added to get the true error. The iterative method applied in 
Section [6T21 was designed to reduce those deviations, but cannot make them zero. 

A better way to determine the real proper motion uncertainties, is to compute the 
standard deviation of the final error- weighted average proper motions. Figure [6] shows these 
scatter-based proper motion errors, for a random sample of about 5% the size of the whole 
catalog. These errors can still be an underestimate of the real proper motion uncertainty, 
because in the error- weighted average the data points are not independent, two overlapping 
bricks reduced into the same plate can produce (positively) correlated proper motions. In 
this case, the true uncertainty is larger than the measured one. 

The most reliable assessment of the proper motion errors is obtained from a comparison 
with external catalogs. The scatter observed in the differences between our proper motions 
and those from Hipparcos, is the combined result of both catalogs' proper motion errors. 
Given that Hipparcos errors are about 1 mas yr -1 , the measured dispersion indicates that 
our real proper motion uncertainties are about 2.3 mas yr -1 , for stars brighter than V = 10. 
This coincides well with the scatter-based errors in Figure EJ at the bright end. 
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The scatter in the proper motion of the LMC and SMC samples can be used to estimate 
the proper motion uncertainties, at the mean magnitude of the Clouds, since their intrinsic 
internal velocity dispersion^ makes a very minor contribution to the observed scatter. Results 
indicate that our real proper motion uncertainties are about 3.8 mas yr _1 for stars around 
V = 16.4, entirely consistent with the scatter-based uncertainties in Figure |6j 

At the faint end, the dispersion in the proper motion of external galaxies varies sub- 
stantially, depending on whether they have CCD or plate 2nd-epoch data. In the first case, 
the dispersion is 11 mas yr _1 at a mean magnitude of V = 17.3. In the second case, the 
dispersion is 21 mas yr _1 at a mean magnitude of V = 18.5. Since the plate photometry 
for the galaxies produced systematically fainter magnitudes than the CCD photometry, the 
difference in these two sets reflects also the difference in precision between (better) CCD 
and (worse) plate measurements, and not just the increase of errors with magnitude. At 
these magnitudes, our scatter-based proper motion uncertainties are about 50% below these 
values. 

In general, the scatter-based proper motion errors are a rather good indicator of the 
real uncertainties in the proper motions for stars brighter than about V=16.5. More im- 
portantly, all the above external estimates: 2.5, 3.8 and 11 mas yr _1 at V ~ 10, 16.4 and 
17.3, respectively, are smaller in size than what was achieved in SPM3, clearly showing the 
increased precision due to having 2nd-epoch CCD data. 



7. Proper Motion of the Magellanic Clouds 

7.1. Selection of LMC and SMC dominated samples 

A photometric selection was made to choose bona fide red giant LMC and SMC stars, 
based on the analysi s of th e 2MASS LMC infrared color magnitude diagram (CMD) of 
Nikolaev &: Weinberg! ( 120001 ) (their sample "J"). The photometric cuts applied, as seen in 
Figure are: 

LMC : 1.1 < J - K s < 1.3 and 9.5 < K s < 12 
SMC : 1.0 < J - K s < 1.2 and 10 < K a < 12.5 

Only stars with CCD 2nd-epoch data were selected, as stars whose proper motions were based 
on plate data only showed a significanly higher dispersion and some visible systematics. 
We were also forced to discard an area of high stellar density, at 71° < a < 76° and 



3 30 km s 1 (0.13 mas yr 1 at 50 kpc distance) or less in each (| Westerlundl 1 1 9971 ) . 
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—68° < 5 < —71°, with CCD data close to the LMC center, which consistently showed 
deviant results probably caused by misidentifications. 

3822 LMC and 964 SMC stars were selected, as seen in Figure El to measure the mean 
absolute proper motion of the Clouds. Bluer sequences of the CMDs in Figure [7] containing 
LMC/SMC stars have signficant contribution from Milky Way stars and therefore were 
not considered. The redder sequences of LMC/SMC AGB stars have very faint V ~ 17.5 
magnitudes, consequently the proper motion errors are too large to be useful. Also, given 
the magnitude-related problems in the plates, it is desirable to have the smallest possible 
difference in brightness between the reference stars and the Clouds' stars, and indeed, the 
chosen samples overlap sufficiently in magnitude (See Figure [9]). 



7.2. Absolute proper motion of the LMC and SMC 



Probability plots (lHamakerlll978l ) of the chosen samples yielded the mean and dispersion 
values for the LMC and SMC proper motion listed in Table [2j The errors quoted include: 
the formal error of the mean value (a j V N stars ) , the error of the quadratic polynomial at 
the LMC and SMC centers, transformation to Hipparcos errors (Eqs. [I] and [2]), and the 
estimated Hipparcos systematic error (0.25 mas yr" 1 ). As explained before, the error budget 
is dominated by Hipparcos systematics. 

Table [3] and Figures [10] and [TTJ summarize how our results compare with recent mea- 
surements of the proper motion of the Magellanic Clouds. The error bars hinder a more 
precise conclusion about the individual tangential velocities of the Clouds based on our 
data. Nonetheless, the methodology used to measure the stellar proper motions in our cata- 
log permits us to make a rather precise measurement of the proper motion of the SMC with 
respect to the LMC, as explained in more detail in Section 17.31 



7.2.1. Center- of-mass proper motion 

The large extent of the Clouds and their non-negligible depth means that all previous 
investigations, which measured proper motions on scattered small fields, had to convert their 
measured values into a center- of-mass proper motion. That is because a given space velocity 
at a fixed distance projects differently on radial velocity and proper motion at different 
locations in the sky, following the same principle of the Moving Cluster method and the 
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Solar Motiorfl. Besides, proper motion obviously scales with dist 



ance. 



Additionally, internal rotation must also be taken into account for the LMC. In the 
case of the SMC, such correction is deeeme d unnecessary because i ts stellar component 



is mostly supported by velocity dispersion ( Harris fc Zaritsky 
tion curve is obtained from radial velocities of Carbon stars ( 



2006). The LMC's rota- 



van der Marel et al. 2002 



Olsen fc Massevl 120071) and yield widely accepted values of V rot) LMC — 50 
Nonetheless, iPiatek et al.l (120081 ) estimates its own V rotj LMC — 120 km s" 1 , 



60 km s 



based on the 

gradient of their measured proper motions along the radius in the LMC disk, and use such 
value for the rotation correction. 



Give n such discrepant v alues, I Cost a et al.l (120091 ) actually obtains two final results for 
the LMC. ICosta et al.l (12009 )-(T) refers to their final proper motion when using V rott LMC — 50 
km s _1 , while ICosta et al.l ( )2009l )-(2) does so for V rot ^MC — 120 km s _1 . For the latter, it 
must be noted that fi Q cos S deviates noticeably from the other determinations in Table |3j 

Altogether, the typical correction for perspective effect for the LMC from both methods 
is about ± 0.2 mas yr" 1 and may rise to about ± 0.5 mas yr" 1 for locations farther than 
5° from the LMC center, running more or less in opposite directions at opposite locations 
on the Cloud. Rotation effects for the LMC are usually less than 0.1 mas yr -1 . Perspective 
effects for SMC are much smaller. For previous studies, these corrections are necessary and 
in some cases yield quite different values for the same fields, but in our work, given the 
spatial extent and symmetry of the data, the net effect on the mean motion of the Clouds is 
very close to zero, and no correction is done. 



7.3. Relative proper motion of the SMC with respect to the LMC 

As explained in Section [672| we measured the mean motion of our reference stars precisely 
all over our field of view, in particular at the location of LMC and SMC within 0.03 and 0.06 
mas yr -1 , respectively. Combined with the relative proper motion of LMC and SMC stars 
with respect to these reference stars, we can indeed measure the proper motion of the SMC 
with respect to that of the LMC, with a higher precision, limited by the error just quoted 
plus the formal error of the mean coming from the number of stars and their measured 
scatter. 

From Tableland taking the errors quoted above into account, it is straightforward to 



7 One of the first attempts to estimate th e proper motion of the Magellanic Clouds was done by measuring 
gradients in the radial velocity along them ( Feitzinger et al. 1977 ; Meatheringham et al. 1988h 
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obtain the relative proper motion of the SMC with respect to the LMC 



Afi acosS (SMC-LMC) 
Afi s (SMC - LMC) 



-0.91 ±0.16 mas yr" 1 
-1.49 ± 0.15 mas yr" 1 



(3) 
(4) 



These values cannot be transformed directly into a measurement of the relative velocity 
between the Clouds. Being at different locations in the sky means that the planes of their 
tangential velocities are different as well, and the necessary rotation and projections to 
measure the SMC velocity on the LMC reference frame does not allow us to obtain the 
relative space velocity as merely a function of the relative proper motion between the Clouds. 

But, we can use these values to obtain new independent measurements of the SMC's 
absolute proper motion, based on existing measurements of the LMC's absolute proper 
motion plus our precise relative proper motion from above. Moreover, since all authors that 
directly measured the proper motion of the SMC had previously measured the LMC's proper 
motion as well, we can verify if their original SMC results are consistent with our relative 
measure. Figure [T2l shows the absolute proper-motion determination for the SMC from Table 
|4] (made by combining our relative LMC-SMC motion with absolute LMC motions from the 
literature), compared to the direct determinations of the absolute proper motion of the SMC 
from Table [3J 



Except for lPiatek et al.1 ( 120081 ) and ICosta et al.l (j2009l )-(2). all published measurements 
of the SMC proper motions are consistent with our new ones. These two works are the 
only ones that use J^l^ mc — 120 km s -1 . It i s wor th noting that field Lll in the LMC 



from iPiatek et al.1 (120081 ) is also in ICosta et al.l (120091 ). but their measured proper motions 
are significantly different, beyond the quote d errors. AH this m akes us suspect that the 
rather small quoted proper motion errors in IPiatek et al.1 (120081 ) underestimate their real 
uncertainties. 

In summary, although our proper motions for the LMC an d t he SMC separately a re in 
agreement - within error bars - with iKallivayalil et al. rt2006bl laL IPiatek et all (boost ) and 



Costa et al] (2009), our relative pr oper m otion of SMC with respect to LMC is consistent 



only with IKallivayalil et~aH (l2006al lbh and lCosta etaJ j2009h -(l). 



7.4. The space motion of the clouds 

The individual as well as relative space velocities of the Clouds, as derived from our 
LMC and SMC proper motions, are given in Table [51 The escape velocity at the dis- 
tance of the LMC is estimated to be 300-350 km s^ 1 , depending on the Galactic potential 
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model used (either a simple isothermal sphere or a more elaborate "cosmologically inspired" 



Navarro- Frank- White (NFW) dark matter profile ([Gardiner &: Noguchilll996t iKlypin et al. 



20021 ; Kue et al.ll2008l ). Taken at face value, the galactocentric velocities in Table indicate 
that the LMC is traveling at a speed that is very close to the escape velocity, while the SMC 
is still below the escape velocity of the Galaxy. Unfortunately, the uncertainties hinder a 
more definitive conclusion regarding their binding status. 

From our proper motion measures, we determine that the SMC is moving at 89 ± 54 
km s _1 with respect to LMC. Our error bars do not allow us to determine whether or not 
the Clouds are bound to each other, given that the escape velocity from the LMC at the 
SMC location is about 90 km s _1 (assuming a simple point-mass geometry, a mass for the 
LMC of 2 x 1O 1O M and 23 kpc for the distance between the Clouds). Yet, we can use 
our new SMC proper motions to obtain additional determinations of the relative velocity 
between the Clouds, using the more precise LMC proper motions available in the literature, 
and the other needed input parameters listed in Table [5j Then we obtain more estimates of 
\\(U,V,W)smc-lmc\\, which are listed in Table [6j For comparison, this table also lists the 
relative velocity originally quoted by the references used. 

It is important to be aware of how sensitive \\(U, V, W)smc-lmc\\ is for different values 
of the LMC proper motion (See Figure [T3|) . It is clear that even the smallest proper motion 
errors quoted translate into a substantial relative velocity error between the Clouds. On top 
of that, the measured values are close enough to the escape velocity of the SMC with respect 
to the LMC, such that any conclusion regarding the binarity of the Clouds is still far from 
decided. The same situation applies to the individual space velocities of the Clouds, the 
binding status of the Clouds to the Milky Way is also extremely sensitive to the individual 
proper motions. 

Thus, although our individual absolute proper motions for the Clouds cannot tell us 
much about their orbits, our measured relative proper motion between LMC and SMC al- 
lowed us to identify the best self-consistent measurements of the proper motion of both 
clouds, from those authors that had measured both galaxies with claimed very good ac- 
curacies. Our measured space velocity has also helped to check which results were more 
consistent with that value, since it is quite sensitive to the proper motions used. 



8. Implications of our results to the dynamics of the Magellanic System 



The most prevalent scenario for the orbit of the Magellanic Clouds, before the HST 
proper motion results, favored them as a binary system in a bound orbit around the Milky 
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Way. Dynamical simulations were performed under the presumption that the Clouds were 
bound to each other, due to the existence of the common HI envelope surrounding them 



( IHindman et al.lll963l ). Crude timing estimates suggested that for its creation and survival, 
a long time of shared orbits was needed. The extent of the Magellanic Stream implied that 
the Clouds had undergone multiple orbits gravitationally bound to the Milky Way and to 
each other. 

The orbits of the Clouds are naturally a key ingredient in such dynamical modeling. 
Once an orbit is chosen, then a full simulation of the Clouds themselves is run. The main goal 
has been to reproduce as much as possible the position and radial velocity of the Magellanic 
Stream, though significant validation is obtained if the Leading Arm, the Magellanic Bridge, 
and the distorted structure of the Clouds, are also replicated. Another "condition" imposed 
is that the Clouds share the same orbital plane with the Magellanic Stream, since the latter 
runs more or less on a great circle on the sky. 

Depending on the particular questions under study, Cloud model investigators use differ- 
ent approaches. Either one or both Clouds are considered, using collisionless (stars) and/or 
collisional (gas) particles. The Cloud(s) is(are) made of one or many particles, which can be 
massless test particles under a given potential, or a conglomerate of self-gravitating particles 
or even "sticky" self-gravitating particles (to model hydrodynamic processes). 

Overall, two competing scenarios have been systematically studied and have had some 
degree of success in the modeling of the mentioned structures: the tidal model and the ram- 
pressure model. In the first scenario, the Galaxy extracted a tidal plume from the LMC 
and/or the SMC, which gave origin to the Magellanic Stream, in a previous close encounter 
of the three bodies about 1.5 Gyr ago. A most recent encounter 200 Myr ago created the 
Magellanic Bridge. The second scenario proposes that the Stream and the Leading Arm 
consist of material that has been ram-pressure-stripped from the LMC (and SMC), during 
its last passage through the extended ionized Halo of the Galaxy, about 500 Myrs ago. 

A complementary scenario to extract subst antial amounts of gas and no stars from the 



Clouds, the blowout model, has been proposed by lNidever et al.l ( 120081 ). The main hypothesis 
is that star formation outflow in the leading edge of the LMC has been blowing out or puffing 
up the gas over the past 2 Gyr, making it easier for ram pressure and tidal forces to strip it 
off. 



The first dynamical simulation of the Magellanic System ( iMurai fc Fujimotolll980l ) was 
done knowing only the radial velocity of the Clouds, and very little about the mass of 
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the Milky Way at the Clouds' distances. Some dynamically permissible parameter^ were 
assumed in order to have all the input parameters necessary to compute an orbit. Their 
models required the inclusion of a massive halo (an idea that was just starting to being 
accepted then) in order to reproduce the highly negative radial velocities observed in the 
tip of the Magellanic Stream. Conversely, assuming that the LMC was in a bound orbit 
provided estimates consistent with a massive Galactic halo. 
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given constant rotational velocity at a large galactocentric distance, while three of the most 
recent calculations use NFW dark matter halo profiles. Not surprisingly, the mass (profile 
and amount) of the Galaxy is an important source of uncertainty in the orbital models. The 
masses of the Clouds are as well another source of error, their current distorted state makes 
any dynamical estimate of the mass a difficult task. 

Given the large mass of the Clouds, dynamical friction^] can significantly reduce their 
perigalacticon distance as they move through the Galactic Halo. N-body simulations of self- 
gravitating particles naturally consider it by d efault, while in othe r case s, it is accounted 



for by using an analytical expression given by iBinnev fc Tremaind ( 119871 ). Recent studies 



( [Hashimoto et al.ll2003t lJust fc Penarrubial 120051 ) have found though, that the latter tends 



to circularize orbits to excess when compared to equivalent N-body simulations, thus some 
simulations have scaled down its effects. 

The important point to consider here is that a high orbital eccentricity, or equivalently 
a high transverse motion with respect to the Galaxy, is needed to lead to the formation of 
the high-velocity Magellanic Stream. Therefore, the initial conditions of the Clouds' orbits 
must be such that even after the effects of dynamical fri ction, orb i tal ec centricity is high 
enough for the Stream to be formed. Interestingly enough, iLin et al.l (jl995[ ) had pointed out 
that a hyperbolic encounter of the Clouds with the Galaxy, in which they are passing by 
for the first (and only) time, could lead to the tidal stripping of gas segments that would 



8 By specifying the orbital inclination of the LMC with respect to the Galactic Plane and its perigalactic 
distance, the equations of energy and angular momentum conservation may be solved to yield the three 
components of the LMC space velocity. 

9 The "retardation" of a moving object when it passes through a region with no n- vanishing mass den sity, 



caused by its gravitational interaction with the particles of that continuous mass (|Chandrasekharlll943l ). 
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later infall and trail rapidly behind. Such a model was discarded then, on the basis that the 
Clouds were presumed bound to the Galaxy. 

Another source of dynamical friction, the LMC dark mat t er ha lo, is considered as well 
in the calculations of the Magellanic orbits ( IBekki fc Chibal 120051 ) . Its effects are quite 
important on the binarity of the galaxy pair, exerting a significant frictional drag on the SMC 
when it penetrates the LMC halo. If they get close enough, they would merge quickly, and 
since this has not occurred yet, it implies either that this force is negl i gible or that the Clouds 
became bound to each other relatively recently ( iGardiner et al.lll994l ). IBekki fc Chibal (120051 ) 
also found that under these conditions, the Clouds cannot keep their binary status for more 
than «5 Gyr in the past. This opens the possibility of the Clouds being coupled only recently. 



Curiously, the very first dynamical model of the Magellanic System (IMurai fc Fujimoto 



19801 ) had to make an extensive search for the binary state that could produce the Magellanic 
Stream. The choice of tangential velocities for SMC had to be so specific, that they could 
- theoretically - predict it within ± 5 km s _1 . This prompted them and others to make the 
first inferred estimates of the proper motion of the Clouds. This early model though, did 
not include a dark matter halo for the LMC. 



A similar situation was faced later on by iHeller fc Rohlfsl ( 119941 ). who explain that if 
the perturbing forces acting on the LMC and/or SMC have to be small enough to leave the 
binary system intact while simultaneously producing long streams similar in shape to the 
Stream, then the evolution time of these streams has to be rather long and very special, and 
properly chosen initial conditions for the test particles in the simulations have to be adopted. 

One point in which all orbital models agree, is that the Clouds had a recent encounter, 
sometime between 200 to 500 Myr ago. In fact, several models also found that the binarity 
of the Clouds, regardless of how long it has been in place, was most probably broken at this 
last collision, which happened very close to their perigalacticton distance. In other words, 
the Milky Way's powerful gravitational tide has disrupted the pair. 

In general, searching for the appropriate binary-bound orbits for the Clouds that can 
later be used to reproduce the Magellanic System, has been a difficult fine tuning task. The 
challenge is now even harder, since the most re cent HST measure ments of the LMC and 
SMC proper motion ( Kallivayalil et al. 2006bl lal IPiatek et al. 2008 ) seem to indicate that 
the Clouds are traveling too fast to have ever been bound to the Milky Way. 



Numerical simulations of the LMC orbit by iBesla et al.l ( 120071 ). based on these new 
numbers, suggest that the Large Cloud is "plunging" in a highly eccentric parabolic orbit, 
on its first passage about the Milky Way. At such speed, dynamical friction is negligible, 
but the choice of the Galactic potential (isothermal sphere vs. NFW) introduces dramatic 
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changes in the orbital history of the LMC. In an isothermal sphere, the LMC has indeed 
a bound orbit, although with an increased period and apogalacticon distance compared to 
previous models. In an NFW profile, the "best case scenario" of a bound orbit has a period 
of about a Hubble time, and reaches an apogalacticon distance of ~ 550 kpc. 

Despite such 3D differences, their projected orbits on the sky were in good agreement, 
so their predicted location compared to the Magellanic Stream's great circle, could not be 
used to distinguish Milky Way mass profile models. More importantly though, in both cases 
the projected orbit did not trace the Magellanic Stream, deviating from it by about 10° on 
the sky. Adding the SMC into the calculations did not reduce the dis agreement, nor using 
the weighted average of all pre-HST proper motions. iBesla et al.l (120071 ) argues then that the 
usual criteria to validate a dynamical model of the Magellanic System, its ability to match 
an orbit with the Stream, is no longer acceptable. 



The conclusion of IBesla et al.l ( 120071 ) , that the Clouds are on their first pas sage relies 
heavil y on the large value for fi a cos 5 that was measured with HST for the LMC. iLin et al. 
( 119951 ) and others have in fact argued that it is hard to explain how a bound LMC could 
have ma ntained a high angular momentum pe r pendicu l ar to the G alaxy's rota t ion ax is for 
so long. iRaychaudhury &: Lynden- Belli ( 119891 ). IShuterl ( 119921 ) and iByrd et al.l (119941 ) have 
suggested that an early interaction with M31 could be the source of such a high tangential 
velocity. In any case, such a condition is easier to understand with an LMC not bound 
to the Milky Way. In ad dition, the Clouds a re the only gas-rich dwarf galaxies at small 
galactrocentric distances ( van den Bergh 20061) . different from the rest in the Local Group. 



All these reasons are used by IBesla et al. 



( 120071 ) to support the case of an unbound LMC. 



Notably, all LMC proper motion results after iKallivayalil et al.l (l2006bl ). including this 
work, have produced lower values of /i Q cos<5. Therefore, a bound orbit with a long period 
and a large apogalacticon distance is still a scenario compatible with the most recent results. 
As for the binarity of the Clouds, a recent period of joint orbits could be enough to explain 
the common features between the Clouds, that is the HI envelope, the Bridge and even the 
Stream, and can account for the different star formation history and chemical evolution in 
each Cloud, that point to a separate origin and place of birth. 

To conclude, the search for a realistic orbit of the Magellanic Clouds is far from over. 
The space velocities obtained in this investigation are supportive of a scenario in which the 
Magellanic Clouds are possibly currently unbound from each other, with the LMC traveling 
at a velocity that is high enough to make it nearly unbound to the Galaxy. But having 
nowadays a much more precise measurement of the proper motion of the Clouds, has not 
facilitated our understanding of their dynamics and has instead opened new questions and 
placed all the constructed models in doubt. 
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9. Conclusions and Future Work 

A catalog of absolute proper motions containing 1,448,438 objects has been obtained 
from SPM material, supplemented with UCAC2 data. The catalog covers an estimated area 
of 450 square degrees except for the inner regions of LMC, SMC and 47 Tucanae, where the 
high stellar density made it impossible to obtain accurate cross-identification of the stars. 

Samples of 3822 LMC stars and 964 stars were selected from the catalog to measure the 
mean proper motion of the Magellanic Clouds. The results obtained are: 

(fia cos 5, us) lmc = (+1.89, +0.39) ± (0.27, 0.27) mas yr" 1 

(fi a cos 5, fi s ) SMC = (+0.98,-1.01) ± (0.30,0.29) mas yr" 1 

Our much more precise relative proper motions with respect to the photometrically 
selected Galactic Disk dwarf stars, enabled us to obtain the proper motion of the SMC with 
respect to the LMC, with significantly smaller uncertainties: 

(/i aC os5, Hs)smc-lmc = (-0.91, -1.49) ± (0.16, 0.15) mas yr _1 

This was used to obtain new independent and more precise proper motions for the SMC, 
based on the more accurate LMC proper motions of other authors. It was also used to 
confirm if their separate measurements of the SMC proper motion were consistent or not 
with our results. 

After a comparison in an absolute and relative sense with previous proper motion results, 
followed by a discussion of the orbital models of the Magellanic Clouds based on those results, 
we conclude that our proper motions are compatible with the LMC and SMC being born 
and formed as separate entities, which later joined in a temporary binary state for the past 
few Gigayears, being recently disrupted by the Milky Way in their most recent perigalaticon 
passage about 200 Myr ago. The Clouds orbits are marginally bound to the Milky Way, 
possibly following a very elongated but still periodic orbit around the Galaxy. 

The search for a realistic orbit of the Magellanic Clouds is far from over. Having 
(formally) very accurate and precise space-based proper motions for the Clouds, has not 
facilitated our understanding of their dynamics but has, instead, opened new questions and 
has placed all dynamic scenarios of the Magellanic System in doubt. As of today, it is 
still unclear if the Magellanic Stream and the Leading Arm are caused mostly by a tidal 
interaction or are the result of the ram-pressure of the Galactic Halo on the gas of the 
Clouds. 

Given the inherent difficulties in measuring an accurate proper motion for the Magellanic 
Clouds, the obvious dangers that systematic errors pose in those measurements and the 
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fact that the dynamical models of the Magellanic System are extremely sensitive to small 
variations in the proper motions of the Clouds, we believe that we are not yet in the position 
of considering them known parameters in the orbital calculation. But we are getting closer. 

9.1. Future work 

From the very begining, this investigation of the proper motion of the Magellanic Clouds 
was known to be restricted by the conditions and characteristics of the SPM material. We 
have proven that even under those constraints, these data are able to produce independent 
significant results on the proper motion of nearby dwarf galaxies. Thus, there is certainly 
room for improvement. 

Second epoch CCD data over the whole field are needed, especially in the inner areas 
of the Clouds. Since the definition of an adequate uncontaminated relative reference frame 
in those areas is difficult, it is also necessary to devise an adequate proper motion reduction 
method, to precisely link these fields into the general global relative reference frame. A 
scaled-down version of the block adjustment method could be put in place here, in which 
plate-size fields of view are assembled first, and then another solution is run to paste those 
into the larger global system. 

Since the dominant limiting factor in terms of precision is the plate measurement errors, 
additional improvement can be achieved by re-scanning the first epoch plates used in this 
work with the Yale PDS, which yields 2 times more precise positions than the currently 
measured ones. As PDS scanning of a full plate is very time-consuming, a subset of properly 
pre-selected objects should be measured. 

Although additional HST follow up observations are already planned to improve the 
space-based proper motions, and other research teams are still working on additional CCD 
ground-based measurements, these studies are still limited since they must correct their 
observed proper motions into a center-of-mass value. Therefore, it is still worth trying to 
improve our results, since they offer a wide field extended coverage of the Magellanic Clouds. 

Another future work being considered, is to search in the Intercloud region for coherent 
structures in the relative proper motion space, to identify stars whose motion is directed 
towards the LMC. This will require to find what possibly is a small number of stars spread 
over an extended area that share systematic proper motions. 

Vieira would like to thank Dr. Burt Jones for being the external reader of her PhD 
thesis, on which this manuscript is based. 



-29 - 



We would like to thank the many former and present members of the Cesco Observatory, 
Universidad Nacional de San Juan and Yale Southern Observatory, who have contributed to 
the SPM program over the years. 

Vieira would also like to thank the National Science Foundation (grants AST04-07292, 
AST04-07293, AST09-08996) and the Yale Astronomy Dept for their financial support during 
her graduate career. She and the other authors are grateful to the NSF and to Yale University 
for support of the SPM and this research in particular throughout the many years required 
to complete the SPM. 

This publication makes use of data products from the Two Micron All Sky Survey, 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/ California Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. 

REFERENCES 

Anguita, C. 1999, in IAU Symposium, Vol. 190, New Views of the Magellanic Clouds, ed. 
Y.-H. Chu, N. Suntzeff, J. Hesser, & D. Bohlender, 475 

Anguita, C, Loyola, P., k Pedreros, M. H. 2000, AJ, 120, 845 

Bastian, U., Roser, S., & Kroupa, P. 1993, in Astronomische Gesellschaft Abstract Series, 
Vol. 8, Astronomische Gesellschaft Abstract Series, ed. G. Klare, 155 

Bekki, K. & Chiba, M. 2005, MNRAS, 356, 680 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

Besla, G., Kallivayalil, N., Hernquist, L., et al. 2007, ApJ, 668, 949 

Binney, J. & Tremaine, S. 1987, Galactic dynamics, ed. S. Binney, J. & Tremaine 

Byrd, G., Valtonen, M., McCall, M., & Innanen, K. 1994, AJ, 107, 2055 

Chandrasekhar, S. 1943, ApJ, 97, 255 

Cioni, M.-R. L., van der Marel, R. P., Loup, C, & Habing, H. J. 2000, A&A, 359, 601 
Connors, T. W., Kawata, D., & Gibson, B. K. 2006, MNRAS, 371, 108 
Costa, E., Mendez, R. A., Pedreros, M. H., et al. 2009, AJ, 137, 4339 



-30 - 

Dehnen, W. & Binney, J. J. 1998, MNRAS, 298, 387 

Drake, A. J., Cook, K. H., Alcock, C, et al. 2001, in Bulletin of the American Astronomical 
Society, Vol. 33, Bulletin of the American Astronomical Society, 1379 

Eichhorn, H. 1960, Astronomische Nachrichten, 285, 233 

Eichhorn, H. 1988, ApJ, 334, 465 

Feitzinger, J. V., Schmidt-Kaler, T., & Isserstedt, J. 1977, A&A, 57, 265 

Freedman, W. L., Madore, B. F., Gibson, B. K., et al. 2001, ApJ, 553, 47 

Gardiner, L. T. & Noguchi, M. 1996, MNRAS, 278, 191 

Gardiner, L. T., Sawa, T., & Fujimoto, M. 1994, MNRAS, 266, 567 

Girard, T. M., Dinescu, D. L, van Altena, W. F., et al. 2004, AJ, 127, 3060 

Girard, T. M., Platais, L, Kozhurina-Platais, V., van Altena, W. F., & Lopez, C. E. 1998, 
AJ, 115, 855 

Hamaker, H. C. 1978, Applied Statistics, 27, 76 

Harris, J. & Zaritsky, D. 2006, AJ, 131, 2514 

Hashimoto, Y., Funato, Y., & Makino, J. 2003, ApJ, 582, 196 

Heller, P. & Rohlfs, K. 1994, A&A, 291, 743 

Hindman, J. V., Kerr, F. J., & McGee, R. X. 1963, Australian Journal of Physics, 16, 570 
H0g, E., Fabricius, C, Makarov, V. V., et al. 2000, A&A, 355, L27 

Irwin, M. J., Demers, S., & Kunkel, W. E. 1996, in Bulletin of the American Astronomical 
Society, Vol. 28, Bulletin of the American Astronomical Society, 932 

Irwin, M. J., Kunkel, W. E., & Demers, S. 1985, Nature, 318, 160 

Jefferys, W. H. 1979, AJ, 84, 1775 

Jones, B. F., Klemola, A. R., & Lin, D. N. C. 1989, in Bulletin of the American Astronomical 
Society, Vol. 21, Bulletin of the American Astronomical Society, 1107 

Jones, B. F., Klemola, A. R., & Lin, D. N. C. 1994, AJ, 107, 1333 



-31 - 



Just, A. & Pefiarrubia, J. 2005, A&A, 431, 861 

Kallivayalil, N., van der Marel, R. P., & Alcock, C. 2006a, ApJ, 652, 1213 

Kallivayalil, N., van der Marel, R. P., Alcock, C, et al. 2005, in Bulletin of the American 
Astronomical Society, Vol. 37, Bulletin of the American Astronomical Society, 1345 

Kallivayalil, N., van der Marel, R. P., Alcock, C, et al. 2006b, ApJ, 638, 772 

Kallivayalil, N., van der Marel, R. P., Anderson, J., Besla, G., & Alcock, C. 2009, in IAU 
Symposium, Vol. 256, IAU Symposium, ed. J. T. van Loon & J. M. Oliveira, 93-98 

Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 

Kroupa, P. & Bastian, U. 1997a, New Astronomy, 2, 77 

Kroupa, P. & Bastian, U. 1997b, in Astronomische Gesellschaft Abstract Series, Vol. 13, 
Astronomische Gesellschaft Abstract Series, ed. R. E. Schielicke, 77 

Kroupa, P. & Bastian, U. 1997c, in ESA Special Publication, Vol. 402, Hipparcos - Venice 
'97, 615-616 

Kroupa, P., Roser, S., & Bastian, U. 1994, MNRAS, 266, 412 
Lee, J.-F. & van Altena, W. 1983, AJ, 88, 1683 

Lin, D. N. C. 1993, in Bulletin of the American Astronomical Society, Vol. 25, Bulletin of 
the American Astronomical Society, 783 

Lin, D. N. C, Jones, B. F., & Klemola, A. R. 1995, ApJ, 439, 652 

Lin, D. N. C. & Lynden-Bell, D. 1982, MNRAS, 198, 707 

Mastropietro, C, Moore, B., Mayer, L., Wadsley, J., & Stadel, J. 2005, MNRAS, 363, 509 

Mathewson, D. S., Cleary, M. N., & Murray, J. D. 1974, ApJ, 190, 291 

Meatheringham, S. J., Dopita, M. A., Ford, H. C, & Webster, B. L. 1988, ApJ, 327, 651 

Mendez, R. A., Costa, E., Pedreros, M. H., & Gallart, C. 2006, in Revista Mexicana de As- 
tronomia y Astrofisica, vol. 27, Vol. 26, Revista Mexicana de Astronomia y Astrofisica 
Conference Series, 183 

Momany, Y. & Zaggia, S. 2005, A&A, 437, 339 



-32 - 



Monet, D. G., Levine, S. E., Canzian, B., et al. 2003, AJ, 125, 984 
Moore, B. & Davis, M. 1994, MNRAS, 270, 209 
Murai, T. k Fujimoto, M. 1980, PASJ, 32, 581 

Nidever, D. L., Majewski, S. R., & Burton, W. B. 2008, ApJ, 679, 432 

Nikolaev, S. & Weinberg, M. D. 2000, ApJ, 542, 804 

Olano, C. A. 2004, A&A, 423, 895 

Olsen, K. A. G. & Massey, P. 2007, ApJ, 656, L61 

Paturel, G., Vauglin, L, Petit, C, et al. 2005, A&A, 430, 751 

Pedreros, M. H. 2006, in Revista Mexicana de Astronomia y Astrofisica, vol. 27, Vol. 25, Re- 
vista Mexicana de Astronomia y Astrofisica Conference Series, ed. C. Abad, A. Bon- 
giovanni, & Y. Guillen, 43-44 

Pedreros, M. H., Anguita, C, & Maza, J. 2002, AJ, 123, 1971 

Pedreros, M. H., Costa, E., & Mendez, R. A. 2006a, AJ, 131, 1461 

Pedreros, M. H., Costa, E., & Mendez, R. A. 2006b, in Revista Mexicana de Astronomia y 
Astrofisica, vol. 27, Vol. 26, Revista Mexicana de Astronomia y Astrofisica Conference 
Series, 78-79 

Penston, M. V. 1982, The Observatory, 102, 174 

Perryman, M. A. C. & ESA, eds. 1997, ESA Special Publication, Vol. 1200, The HIPPARCOS 
and TYCHO catalogues. Astrometric and photometric star catalogues derived from 
the ESA HIPPARCOS Space Astrometry Mission 

Piatek, S., Pryor, C, & Olszewski, E. W. 2008, AJ, 135, 1024 

Platais, L, Girard, T. M., Kozhurina-Platais, V., et al. 1998, AJ, 116, 2556 

Putman, M. E., Gibson, B. K., Staveley-Smith, L., et al. 1998, Nature, 394, 752 

Raychaudhury, S. & Lynden-Bell, D. 1989, MNRAS, 240, 195 

Shuter, W. L. H. 1992, ApJ, 386, 101 

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 



-33 - 

Stanimirovic, S., Staveley-Smith, L., & Jones, P. A. 2004, ApJ, 604, 176 
Stock, J. 1981, Revista Mexicana de Astronomia y Astrofisica, vol. 6, 6, 115 
Taff, L. G. 1988, AJ, 96, 409 

Tucholke, H. & Hiesgen, M. 1991, in IAU Symposium, Vol. 148, The Magellanic Clouds, ed. 
R. Haynes & D. Milne, 491 

van Altena, W. F., Girard, T. M., Platais, I., et al. 1999, in Bulletin of the American 
Astronomical Society, Vol. 31, Bulletin of the American Astronomical Society, 1229 

van den Bergh, S. 2006, AJ, 132, 1571 

van der Marel, R. P., Alves, D. R., Hardy, E., & Suntzeff, N. B. 2002, AJ, 124, 2639 

van Leeuwen, F., ed. 2007, Astrophysics and Space Science Library, Vol. 350, Hipparcos, the 
New Reduction of the Raw Data 

Veron-Cetty, M.-P. & Veron, P. 2006, A&A, 455, 773 

Wesselink, A. J. 1974, in IAU Symposium, Vol. 61, IAU Symp. 61: New Problems in As- 
trometry, ed. W. Gliese, C. A. Murray, & R. H. Tucker, 201 

Westerlund, B. E. 1997, Cambridge Astrophysics Series, 29 

Wilson, R. E. 1915, Proceedings of the National Academy of Science, 1, 183 

Xue, X. X., Rix, H. W., Zhao, G., et al. 2008, ApJ, 684, 1143 

Yoshizawa, A. M. & Noguchi, M. 2003, MNRAS, 339, 1135 

Zacharias, N., Urban, S. E., Zacharias, M. I., et al. 2004, AJ, 127, 3043 

Zacharias, N., Winter, L., Holdenried, E. R., et al. 2008, PASP, 120, 644 



This preprint was prepared with the A AS IATgX macros v5.2. 



-34- 




Fig. 1. — SPM fields in the Magellanic Cloud area. Gray outlined regions correspond to 
2 nd -epoch SPM plate observations. The small squares indicate 2 nd -epoch CCD observations. 
The black outlined area is the VMC region (See Sec. 14. 2D . A subset of stars from the USNO- 
A2.0 Catalog is plotted in the background, to indicate the positions of the LMC (on the left) 
and SMC (on the right). SPM fields are on nominal 5° centers, not including the VMC field. 
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Fig. 2. — Magnitude distribution and, thus, completeness of the SPM material evaluated in 
a CCD-size field in SPM field 028 that has SPM plate data (dark grey), SPM CCD data 
(light grey) and UCAC2 data (stitch under curve), in comparison with the 2MASS detections 
(black) in the field. 
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Fig. 3. — Sky map of the relative proper-motion reference stars. The grey curves enclose the 
so called contaminated areas. 
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Fig. 4. — Systematics when comparing a single frame with the average frame. This example 
illustrates the case of a PV frame, for which A/ig = H s single frame - average frame exhibits 
visible trends, as shown by the median line. 
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Fig. 5. — Sky coverage of the final catalog of proper motions. Every tenth entry is plotted 
in light gray. The "contaminated" areas are as in Figure EJ The lower density area north of 
the SMC corresponds to an area where several poor quality plates were discarded. 
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Fig. 6. — Uncertainty in [x$ vs. magnitude. The error estimate is based on the scatter of the 
individual data points that were error-weighted-averaged into the final proper motion. The 
black line is a moving median, computed over a 0.5 mag interval. The gray dashed line is 
the moving median value for /i a cos 5. This plot is based on a random sample of 5% the size 
of the whole catalog. 
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Fig. 7.— Infrared CMDs in the vicinity of the LMC (left) and SMC (right) based on 2MASS 
photometry of the stars in our proper motion catalog. The two bluest vertical sequences 
visible in both CMDs are dominated by foreground Galactic stars, while the bright red 
sequences are dominated by LMC/SMC stars. Black points on each planel show the red 
giant star samples selected to measure the mean motion of each Cloud. 




Fig. 8. — Distribution on the sky of LMC and SMC sample stars selected from our catalog 
of proper motions. 
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Fig. 9. — Histogram of V magnitude for the reference stars, LMC and SMC stars. All three 
data sets overlap substantially in V, which helps to reduce systematics in the proper motions 
related to brightness. 
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Fig. 10. — Proper M otion of the Large Magellanic Cloud. Labels are as follows : 



J94= l Jones et all fll994J). K9 7= lKrouDa fc Bastianl (11997dl . K0 6=lKallivavalil et all (12006bh . 
P0 6=IPedreros et al.N2006~l . P08= iPiatek et all fhoosh . C09-l= ICosta et al.l hoQW (1). C09- 
2= ICosta et all (120091 ) (2). All these results are plotted in gray while ours is plotted in black. 
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Fig. 11. — Proper Motion of th e Smal l Magellanic Cloud. Labels are as fol 
l ows: K97= lKroupa & BastiarJ Jl997ch . K06= lKallivavalil et all fj2006ah . P08= |Piatek et al. 
(120081 ) .C09= ICosta et ajj (120091 ) . All these results are plotted in gray while ours is plotted in 
black. 



-43- 




Fig. 12. — Comparison of previous and new determinations of the Proper Motion of the SMC. 
Previous determinations from Figure (TT] are plotted with filled translucent symbols: K97 in 
yellow, K06 in blue, P08 in green, C09 in red. Our new determinations based on the same 
studies' LMC proper motion plus our relative SMC-LMC proper motion, are plotted with 
open symbols of the same colors (C09 has two determinations of the LMC proper motion). 
Previous results are consistent with ours when the filled and the open symbol of the same 
color overlap. Our SMC proper motion in plotted with the black symbol. 
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Fig. 13. — Relative velocity between the Clouds as a function of the absolute proper motion of 
the LMC. For any given value (fi a cos 8, /j,s) lmc hi the diagram, (fi a cos 5, Hs)smc is obtained 
from Eqs. |3] and H] and then \\(U,V,W)smc-lmc\\ is computed (other needed parameters 
from Table EJ) and plotted with a grey color indicative of its value. LMC proper motions and 
errors from K06, P08, C09-1 and C09-2 are plotted in black with their error ellipses, and 
from this work with error bars only. The black parabollic curve indicates the values of LMC 
proper motion for which \\(U, V, W)smc-lmc\\ — 90 km s _1 . Values of LMC proper motion 
that fall below this curve roughly correspond to the Clouds forming a bound binary system. 
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Table 1. SPM plates used in this paper 



SPM plates 

SPM field (a, S) 1950 1st epoch 2nd epoch 









Plate number 


Date 


Plate number 


Date 


028 


(00 h 00 m , 


-75°) 


0829B/Y 


03/NOV/72.84 


1354/BY 


01/NOV/94.84 


029 


(01 h 00 m , 


-75°) 


0834B/Y 


06/NOV/72.85 






030 


(02 h 00 m , 


-75°) 


0826B/Y 


01/NOV/72.84 






031 


(03 h 00 m , 


-75°) 


0830B/Y 


03/NOV/72.84 






032 


(04 h 00 m , 


-75°) 


0950B/Y 


26/OCT/73.82 






033 


(05 h 00 m , 


-75°) 


0851B/Y 


06/JAN/73.02 


1357/BY 
1373/BY 


30/NOV/94.92 
20/OCT/95.80 


034 


(06 h 00 m , 


-75°) 


0860B/Y 


25/JAN/73.07 






052 


(00 h 00 m , 


-70°) 


0821B/Y 


06/OCT/72.77 


1355/BY 


02/NOV/94.84 


053 


(00 h 48 m , 


-70°) 


0750B/Y 


24/AUG/71.65 


1371/BY 


23/SEP/95.73 


054 


(01 h 36 m , 


-70°) 


0751B/Y 


25/AUG/71.65 


1372/BY 


20/OCT/95.80 


055 


(02 h 24 m , 


-70°) 


0969B/Y 


26/NOV/73.90 






056 


(03 h 12 m , 


-70°) 


0815B/Y 


16/SEP/72.71 






057 


(04 h 00 m , 


-70°) 


0946B/Y 


06/OCT/73.76 






058 


(04 h 48 m , 


-70°) 


0846B/Y 


01/JAN/73.00 






059 


(05 h 36 m , 


-70°) 


0854B/Y 


09/JAN/73.02 


1350/BY 


17/NOV/93.88 


060 


(06 h 24 m , 


-70°) 


0842B/Y 


13/DEC/72.95 






061 


(07 h 12 m , 


-70°) 


0423B/Y 


20/FEB/69.14 






088 


(04 h 00 m , 
(04 h 40 m , 


-65°) 


0824B/Y 


15/OCT/72.79 






089 


-65°) 


0835B/Y 


06/NOV/72.85 






090 


(05 h 20 m , 


-65°) 


0853B/Y 


07/JAN/73.02 


1394/BY 


16/NOV/93.88 


091 


(06 h 00 m , 


-65°) 


0781B/Y 
0984B/Y 


13/JAN/72.04 
27/JAN/74.07 






092 


(06 h 40 m , 


-65°) 


0986B/Y 


28/JAN/74.08 







Note. — Plates 0750B, 0751B, 1371B, 1357Y, and 1373Y, were subsequently discarded, as visual 
examination of them revealed significant deficiencies that would negatively impact the astrometric 
reduction at their locations. 
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Table 2. Absolute Proper Motion of the Magellanic Clouds 



Target 




fi a cos 5 
mas yr _1 


US 

mas yr _1 


£fi a cos 5 

mas yr _1 


e V5 

mas yr _1 


COS (5 

mas yr _1 


mas yr _1 


LMC 
SMC 


3822 
964 


1.89 
0.98 


0.39 
-1.01 


0.27 
0.30 


0.27 
0.29 


3.76 
4.12 


3.59 
3.82 



Note. — The e values are the errors of the mean \i values, which include: the formal 
errors (cr/N s tars), the error of the quadratic polynomial at the LMC and SCM centers, 
transformation to Hipparcos errors (Eqs. [T]and|2j), and the Hipparcos systematic error 
(0.25 mas yr _1 ). The a values represent the scatter of the data around the mean proper 
motion values. 
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Table 3. Recent Determinations of the Proper Motion of the Magellanic Clouds 



Author 


fi a cos 5 


fJ-5 


Nstars 




mas yr _1 


mas yr _1 




Large 


Magellanic Cloud 




Jones et al. ( 1994) 


1.37 ±0.18 


0.28 ±0.27 


251 


Kroupa h Bastian ( 1997c) 


1.94 ±0.29 


-0.14 ±0.36 


33 


Pedreros et al. (2006a) 


1.80 ±0.10 


0.90 ±0.10 


108 


Kallivavalil et al. (2006b) 


2.03 ±0.08 


0.44 ± 0.05 


810 


Piatek et al. (2008) 


1.95 ±0.04 


0.44 ± 0.04 


889 


Costa et al. (2009) (1) 


1.82 ±0.13 


0.39 ±0.15 


41 


Costa et al. (2009) (2) 


1.61 ±0.13 


0.60 ±0.15 


41 


This work 


1.89 ±0.27 


0.39 ±0.27 


3822 


Small 


Magellanic Cloud 




KrouDa &; Bastian (1997c) 


1.23 ±0.84 


-1.21 ±0.75 


9 


Kallivavalil et al. (2006a) 


1.16 ±0.18 


-1.17 ±0.18 


177 


Piatek et al. (2008) 


0.75 ±0.06 


-1.25 ±0.06 


215 


Costa et al. (2009) 


1.03 ±0.29 


-1.09 ±0.18 


44 


This work 


0.98 ±0.30 


-1.10 ±0.29 


964 



Note. — Results (1) and (2) from lCosta et al.l (120091 ) are obtained 
assuming two different rotational velocities for LMC, as explained in 
detail in Section ITT3"1 
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Table 4. New determinations of the Proper Motion of the SMC 



Using the LMC 


/J, a COS 5 




proper motion of 


mas yr _1 


mas yr _1 


Jones et al. (1994) 


0.46 ±0.31 


-1.21 ±0.31 


KrouDa & Bastian fl997c) 


1.03 ±0.39 


-1.63 ±0.39 


Pedreros et al. (2006a) 


0.89 ±0.18 


-0.59 ±0.19 


Kallivavalil et al. (2006b) 


1.12 ±0.16 


-1.05 ±0.18 


Piatek et al. (2008) 


1.04 ±0.16 


-1.05 ±0.16 


Costa et al. (2009) (1) 


0.90 ±0.21 


-1.10 ±0.21 


Costa et al. (2009) (2) 


0.70 ±0.21 


-0.60 ±0.21 



Note. — Results (1) and (2) from lCosta et al.1 $200% are 
obtained assuming two different rotational velocities for 
the LMC, as explained in detail in Section 17.31 Compare 
with the direct determinations of the proper motion of the 
SMC in Table Quoted errors include the contribution 
from our measured relative proper motion. 



Table 5. Space velocity of the Magellanic Clouds 



Parameter 


Unit 


LMC 


SMC 




Measured or 


Adopted Quantities 




(a, 8) 


deg 


(81.90,-69.87) 


(13.20,-72.50) 


(l,b) 


deg 


(280.526,-32.527) 


(303.788,-44.628) 


Distance 


kpc 


50.1 ± 2.3 


62.8 ± 2.6 


Radial velocity 


km s _1 


262.1 ± 3.4 


146.0 ± 0.6 


(Ha cos 5, us) 


mas yr _1 


(1.89,0.39) 


(0.98,-1.10) 






± (0.27,0.27) 


± (0.30,0.29) 


(U,V,W) Q , LSR 


km s _1 


(10.00,5.25,7.17) 


(U,V,W) LSR , gc 


km s _1 


(0,220,0) 




Derived Quantities 




(X,Y,Z) 


kpc 


(0.3,-41.5,-26.9) 


(-16.2,-37.6,-44.1) 


Rgc 


kpc 


41.5 


40.9 


(U,V,W) hc 


km s _1 


(71.8,-468.1,226.1) 


(51.5,-433.4,139.2) 






± (63.6,37.1,56.6) 


± (82.1,70.3,62.2) 


(U,V,W) gc 


km s _1 


(71.8,-248.1,226.1) 


(51.5,-213.4,139.2) 






± (63.6,37.1,56.6) 


± (82.1,70.3,62.2) 


\\(U,V,W) gc \\ 


km s -1 


343.3 ± 47.8 


259.9 ± 68.6 



\\{U,V,W)smc-lmc\\ km s- 1 89.2 ± 53.9 



Note. 



Stanimirovic et al 



( a, 6) from Ivan der Marel et all ()200a ) (LMC) and 



(2004) (SMC). (l,b) ar e transforme d from (a, 5). Distance 



from iFreedman et al.1 (1200111 (LMC ) and Cioni et al.1 (1200m (SMC). Radial 
veloci ty from van der Marel et al. ( 20021 ) (LMC) and Harris Sz Zaritskv 



(|2006l ) (SMC). Solar motion (U, V, W) &j lsr, and galactocentric vel o city o f 
the Local Standard of Rest (U, V, W)LSR,gc, from iDehnen &; Binneyl (jl998l ). 
Rq = 8 kpc was adopted. R gc is the Galactocentric distance. (U, V, W)h c 
and (U, V, W) gc are the heliocentric and galactocentric space velocities of the 
Clouds. 1 1 (U, V, W)lmc,SMC 1 1 is the relative velocity of SMC with respect to 
LMC. 
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Table 6. Relative velocity between LMC and SMC 



Reference 


\\(U,V,W)SMC-LMC\\ 




km s -1 


Originally quoted by reference 


Kallivavalil et al. (2006a) 


105 ± 42 


Piatek et al. (2008) 


142 ± 19 


Costa et al. (2009 Wll 


84 ± 50 


Costa et al. (2009W2) 


62 ± 63 


New value using 


our HSMC-LMC 


Kallivavalil et al. (2006a) 


123 ± 51 


Piatek et al. (20081 


108 ± 51 


Costa et al. (2009Vm 


83 ± 52 


Costa et al. (2009V (2) 


152 ± 63 


This work 


This work 


89 ± 54 



Note. 



Results (1) and (2) from ()Costa et al 



20091 ) are obtained assuming two different rotational 



velocities for LMC, as explained in detail in Section 
1731 



